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Abstract 

The year 2000 marks the centenary of the Planck function for the 
blackbody spectral distribution. Here are traced the steps, over 100 years, 
along the path leading to this important formula. 

1 Introduction 

The Planck function was first presented at a meeting of the German Physical 
Society in Berlin on 19 October 1900 in a contribution entitled "On an im- 
provement of the Wien spectral equation" , and published in the proceedings of 
the society in a brief three-page report (Planck, 1900Va). Only the functional 
form was given, including two undetermined constants, derived on the basis of a 
tentative relation between the entropy and the energy of each of a collection of 
radiation resonators. The definitive formula was presented at a meeting of the 
society on 14 December 1900 in a paper entitled "On the theory of the law of 
the energy distribution in the normal spectrum" and published, in a nine-page 
report (Planck, 1900Vb). Herein Planck introduced his new natural constant 
h relating the "energy element" e associated with a resonator of frequency v 
by the formula e = hv, with the numerical value h = 6.55xl0 -27 erg seconds. 
He only sketches the actual calculation of the entropy, now for the first time 
related to the probability of resonator complexions via the Boltzmann princi- 
ple. Thus ended the long search for the blackbody spectral distribution function 
and simultaneously began the era of quantum physics. A full treatment of the 
derivation was later given in a paper (Planck, 1901a) entitled "On the law of the 
energy distribution in the normal spectrum" submitted to Annalen der Physik 
on 9 January 1901. 
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2 Thermal Radiation 



That the arrival of both the visible light and the 'dark' radiant heat from the 
Sun are simultaneously interrupted and restored at the time of a solar eclipse, 
indicated that light and heat are propagated across space with the same veloc- 
ity. The astronomer W. Herschel (1800), using a thermometer exposed to the 
prismatic spectrum of the Sun, discovered that the heating properties of light 
extended beyond the visible spectrum, being enhanced most especially in the 
infrared. He also showed that the unseen heat radiation obeyed the same laws 
of reflection and refraction as visible light. The similarity in other properties, 
including bi-refringcncc and polarisation, as found by others, for light and heat, 
illustrated a difference of degree rather than kind. Furthermore, that sources 
of light are also sources of heat, and vice versa, served to distinguish the two 
manifestations only by the means of detection employed. Indeed, from his mea- 
surements, Herschel presented two distributions of the (continuous) intensity of 
the solar rays, based respectively on their ability to illuminate and to heat. Be- 
cause of the two different techniques of measurement, the relative scales are not 
related. If the two intensity distributions were brought to a common scale and 
then summed, it would represent the first wideband total intensity distribution 
for the solar spectrum. 

Simple experiments, such as those by Leslie (1804) and others, showed that 
the rate of emission of thermal radiation from a heated body, expressed as the 
emissive power or radiant emittance E, being the power emitted per unit area, 
depends on the temperature and on the nature of the surface. Thus a blackened 
surface has a higher emittance than a polished one at the same temperature. 
Also when two such surfaces are subjected to incident radiant energy the fraction 
of energy absorbed, or absorptivity A, depends on the temperature and on the 
nature of the surface, a blackened surface being a more efficient absorber than 
a polished one. 

Kirchhoff (1859), by means of a thought experiment, demonstrated that the 
ratio E/A is the same for all bodies at the same temperature, emitting and 
absorbing at the same wavelength. He considered two bodies, in the form of 
infinte surfaces, Si and S2, backed with perfectly reflecting mirrors and arranged 
parallel to each other with the mirrors on the outside. Si is such that it can 
only emit and absorb radiation of wavelength A, while S2 can emit and absorb 
all possible wavelengths A. At wavelength A the emissivities are E\ and E 2l and 
the absorptivities are A\ and A 2 . For all A ^ A all the radiation emitted by 
S2 is reabsorbed by it. Therefore for temperature equilibrium S2 must absorb 
as much as it emits of the radiation of wavelength A = A, whether it comes 
from Si or S2. Of the radiation E\ emitted by Si, S2 absorbs, after repeated 
absorptions and emissions : 

00 

A 2 E 1 J2 kn = A 2 E 1 /(l-k) 

n=0 



2 



while of the radiation E 2 emitted by S2 



itself, S2 absorbs: 



00 



(1 - A 1 )A 2 E 2 ]T k 



(l-A 1 )A 2 E 2 /(l-k) 



n=0 

where k = (1 - Ai)(l - A 2 ). Hence : 
E 2 = A 2 [E\ + (1 



A l )E 2 ]/{\ - k) 



giving : 



A X E 2 = A 2 E X or 



E 1 /A 1 = E 2 /A 2 



The same result is obtained on applying the condition of temperature equilib- 
rium to Si, yielding : 



Since the composition of the surfaces, as well as the wavelength A and the 
temperature are arbitrary it follows that the ratio E/A, at any wavelength or 
temperature, is the same for all surfaces, and is a function only of wavelength 
and temperature, which statement is known as Kirchhoff's Law. 

Kirchhoff (1860) generalised the above argument, which involved only dis- 
crete wavelengths, to the continuous spectrum, by re-defining E as the emittance 
per unit wavelength interval. He also introduced the concept of a perfectly black 
body, being one which absorbs all incident radiation, for which A = 1 at all 
wavelengths and temperatures. By considering the temperature equilibrium of 
a black body enclosed in a cavity with walls made of a similarly black mate- 
rial, he showed that the emittance is equal to the radiant power, F(A,T) per 
unit wavelength interval, crossing unit area, within the 'black body' enclosure. 
Hence, by considering the temperature equilibrium of a non-black body within 
the cavity, Kirchhoff's Law takes the form E/A = F(A,T). Thus the absorp- 
tivity of a body is equal to its relative emittance, or the emittance of the body 
relative to that of a black body. Kirchhoff stressed the importance of the study 
of black body radiation in order to determine the function F(X,T). It is also 
worth noting here that the enunciation by Kirchhoff of the principle of detailed 
balance between emission and absorption precedes that of Einstein by more 
than half a century. 

3 Integrated Properties of Black Body Radia- 
tion 

Given the dependence on wavelength A and temperature T of the black body 
radiation intensity /(A, T), its energy density u and flux F, both monochromatic 
and integrated, can be obtained : 



E 1 =A 1 [E 2 + {l-A 2 )E l }/{l-k) 
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u(X,T) = (4ir/c)I{X,T),u(T) = J u(X,T)d\ = {4n/c)I(T) 

F(X, T) = nI(X, T), F(T) = J F(X, T)dX = (c/4)«(T) 

As indicated by Kirchhoff, the properties of black body radiation can be 
studied by observations on the radiation emerging from a small orifice in the 
wall of a black body cavity. Experiments by Stefan (1879) established that the 
total (integrated) emergent flux could be represented in the form : 

F(T) = aT 4 

where a — 5.Q70xlO~ s Jm~ 2 s^ 1 K ~ 4 is the Stefan(-Boltzmann) constant, from 
which it follows that the energy density is given by : 

u{T) = (Aa/c)T 4 = aT 4 

with the radiation constant a = 7.5Q5xl0~ 16 Jm~ 3 K~ 4 . The latter form was 
established theoretically by Boltzmann (1884) on the basis of thermodynamics. 
Using for the pressure of radiation the relation : 

P(T) = (1/3)«(T) 

as indicated by the electromagnetic theory of Maxwell (1873), and assuming that 
radiation obeyed the known laws of thermodynamics, then putting U = uV : 

dQ = TdS = d(uV) + PdV = {u + P)dV + Vdu 

from which it follows that : 

(d/du) v [(u + P)/T] = {d/dV) u [V/T] 

which, on using P = u/3, gives : 

(T/u)(du/dT) v = (dlnw/dlnTV = 4 

so that : 

u = aT 4 , P=(l/3)aT 4 

in agreement with the result of Stefan. Alternatively, Stefan's result may be 
taken as a proof of the relation P = u/3. From the above expressions for u and 
P it follows that : 

dS = (4:/3)aT 3 dV + 4aT 2 VdT 

giving for the entropy : 

S - (A/3)aT 3 V 
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4 Spectrum of Black Body Radiation 



The challenge, first presented by Kirchhoff (1859), to determine the spectral 
distribution of black body radiation, was not to be completely answered until 
the turn of the century. Experiments by Draper (1847) and by Weber (1887) 
showed that, with increasing temperature, the maximum intensity in the light 
emitted by incandescent solid bodies, progressed across the spectrum, in the 
direction from the red to the violet. Bartoli (1876) had pointed out that, in 
order not to contradict the second law of thermodynamics, radiation must exert 
a pressure, and Boltzmann (1884) showed that work must be done against this 
pressure in compressing a volume occupied by radiation alone. Wien (1893), 
developing this argument, made the first important steps in the determination 
of the black body radiation spectrum. On the basis of the Doppler principle, 
applied to radiation within a rectangular parallelepiped and reflected (normally) 
from a wall moving with velocity v, where the relative wavelength change is given 
by 5X/ A = 2v/c (twice the ordinary Doppler shift), then radiation of wavelength 
Ao and frequency vq in a volume Vo, subjected to an adiabatic compression to 
volume V, undergoes a change to wavelength A and frequency v where : 

\ 3 V~ 1 = AqV^ 1 and v 3 V = v$V 

The change in energy of the radiation, supplied by the work done in the com- 
pression against the radiation pressure, results in a change in temperature from 
T to T. Since for an adiabatic change VT 3 — VqT^ it follows that : 

AT = A T and vjT = u /T 

giving the earliest expression of the Wien Displacement Law, which states that 
with changing temperature the black body spectrum is shifted so that for each 
wavelength the product of the temperature and the wavelength remains con- 
stant. Similarly each interval of wavelength and frequency is changed according 
to : 

TdA = T dA and du/T = du /T 
Since the (spatial) energy density obeys the relation : 

M (A,T)dA/ U (Ao,To)dA - u(u, T)du/u(u , T )d^ - T 4 /T 4 

then : 

u(A,T)/«(Ao,T ) = (T/T ) 5 = (A /A) 5 

u{v,T)/u{v ,T ) = (T/T ) 3 = {y/y Q f 

in agreement with the observations of Weber. It also follows that the wavelength 
A m (T) and frequency v m {T) at which the wavelength and frequency spectral 
distributions reach their maxima satisfy the relations : 

T\ m (T) = T A m (T ) and u m {T)/T = u m (T )/T 
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representing another expression of the Wicn Displacement Law. Therefore given 
the spectral distribution for one temperature To, that for any other temperature 
T can be determined. 

Wien (1894) re-derived all the above relations when the radiation is enclosed 
in a spherical volume of variable radius. He also showed that the spectral 
distributions for different temperatures cannot cross. The condition for this to 
be true is that : 

[ain«(A,T)/ainA] r > -5 

which is trivially true for the rising branch of the curve, and also obviously 
true for the falling branch. In terms of the frequency distribution the condition 
becomes : 

[d\nu{v,T)/d\nv] T < +3 
4.1 The Wien Formula 

Wien (1896) pointed out that, since the black body radiation spectrum corre- 
sponded to thermal equilibrium, it should be completely determined from the 
principle of maximum entropy, although he did not do this himself. Instead, 
adopting a radiation model based on analogy with the Maxwell velocity distri- 
bution, he obtained a spectrum of the form : 

U (A,T) = / 1 (A)cxp[-/ 2 (A)/T] 

where /i(A) = ci/A Q , (ci and a being constants), and /2(A) are functions to 
be determined. In order to satisfy the Stefan-Boltzmann T 4 dependence of the 
energy density, it is necessary for a = 5 and for / 2 (A)/T to be a function of 
the product AT, thus giving a more general mathematical expression to the 
Wien Displacement Law. Putting, specifically, /2(A)/T = c 2 /(AT), where c 2 is 
a constant, gives the form : 

u(X, T) = ciA~ 5 cxp[-c 2 /(AT)] 

known as Wien's formula, which attains its maximum value at a wavelength X rn 
with A m T = C2/5. The formula was in good accord with the observations of 
Paschen (1896), with c 2 = 5A m T = 14455/U K, or A m T = 2.891xl(T 3 m K. 

In a series of studies, from 1897 onwards, on the entropy of radiation, treating 
the radiation field as an ensemble of monochromatic resonators, Planck sought 
to realize the hope of Wien that the blackbody radiation spectrum could be 
derived on the basis of the maximization of the entropy. Indeed Planck (1899) 
believed (erroneously) that he had shown that the Wien spectral distribution 
was a necessary consequence of the principle of entropy maximization, and that 
its range of validity coincided with that of the (two) laws of thermodynamics. 
Such support for the Wien formula led to it being referred to as the Wien-Planck 
distribution. However, experimental measurements by Lummcr and Pringsheim 
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(1899), Rubens (1899) and others pointed to systematic departures of the Wien- 
Planck formula from the experimental results particularly at long wavelengths. 
Thicsen (1900), following the lead of Wien, expressed the spectral distribution 
in the, only slightly different, general form : 

u(X,T) =T 5 V(AT) 

where tp(XT) is a still unknown function. Thicsen also suggested that the func- 
tion tp, when determined, would accommodate the deviations between theory 
and experiment by differing from the Wien-Planck form, leading to a modifica- 
tion of the Wien displacement law. 

Planck advanced his support of the Wien formula in two further papers 
(Planck, 1900Aa and 1900Ab). In the first paper he adopted, by way of defi- 
nition, s = — (w/av) h\{w/ebv), where s and w are the entropy and (average) 
energy per radiation resonator of frequency v, and a and b are universal posi- 
tive constants, while e is the base of natural (Napierian) logarithms. Using the 
thermodynamic relation ds/dw = 1/T where T is the temperature, it follows 
that w = bv cxp (—av/T), from which, on multiplying by the density Sirv 2 /c 3 
of resonators in ^-space, the Wien formula follows. From the measurements of 
Paschen (1899) the constants a and b were determined to be a = 4.818a;10 _11 
sec K and b = 6.885x10—27 erg sec. The value of b here therefore represents 
the first determination of what later became the Planck constant h, while the 
value of a is that of h/k. In the second paper Planck sought to justify his ear- 
lier assumed definition of the entropy s, by arguing that maximization required 
d 2 s/dw 2 = —F(w), where F(w) is a positive function of w. Using a further rela- 
tion (of limited validity) F(nw) = F(w)/n, where n is the number of resonators, 
Planck adopted the solution F(w) cx 1/w. Thence putting d 2 s/dw 2 = —a/w, 
it follows that s = —awln(@w/e), where now a and (3 are positive constants 
which depend only on the frequency v. Then w = [l/e(3(v)] exp[— 1/Ta(v)], 
leading again to a spectral distribution of the Wien form. Furthermore the 
Wien displacement law required both a and (3 to be proportional to v~ x , thus 
permitting the identifications a = 1/av and [3 — l/ebv, where a and b are 
the quantities introduced earlier. However further work by Lummer and Pring- 
sheim (1900) showed that both a and b were not constant but increased with 
wavelength, leading them to conclude the invalidity of the Wien-Planck spectral 
distribution. 

4.2 The Rayleigh Formula 

To Rayleigh (1900) the Wien formula was difficult to accept on the grounds 
it implied that, with increasing temperature, the radiation density at a given 
wavelength approached a limit. Using an argument based on the still disputed 
Maxwell-Boltzmann principle of the equipartition of energy, in which each de- 
gree of freedom has an energy proportional to the temperature T, and in analogy 
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with the theory of sound, taking for the density of vibrational modes n(X) oc A 4 , 
he obtained : 

u(A,T) =KT\- A 

where K is a constant, which has the form prescribed by Thiesen, suggesting 
that it may be correct for large values of AT. However the formula tends to 
infinity at short wavelengths, leading Rayleigh to attach the Wien exponential 
factor exp[— c 2 /(AT)]. 

4.3 The Planck Formula 

In the light of the increasing evidence for the inadequacy of the Wien blackbody 
spectral distribution, Planck (1900Va) now addressed the problem of improving 
upon it. He recalled his earlier conclusion that the maximization of entropy 
was not of itself sufficient to derive the correct distribution, and recognised 
that his support of the Wien function had depended upon his adoption of the 
further special condition d 2 s/dw 2 = —a/w, applicable only in the case of indef- 
initely small values of w and s. Accordingly he considered instead the condition 
cPs/dw 2 = —a/w(/3 + w), which he described as by far the next simplest rela- 
tion giving s as a logarithmic function of w, while at the same time reducing 
to the former relation for vanishingly small values of w. It is readily shown 
that, ignoring constants of integration, ds/dw = — {a/ (3) ln{w/(/3 + w)}, and 
s = — (a/ (3){w\nw — ((3 + w) ln(/3 + w)}. Then, using ds/dw = 1/T, it follows 
that w = /3/{exp(/3/aT) — 1}. Noting that the Wien Displacement Law requires 
(see later) that s be a function of wjv, where v is the frequency, (demanding 
that (3 be proportional to v), then enabled Planck to express the spectral dis- 
tribution in the form u(X,T) — G\X~ 5 /{exp(C2/ AT) — 1}, where C\ and Ci are 
constants, thus introducing what is now known as the Planck Function. Rubens 
and Kurlbaum (1900) pointed out that the new Planck formula included those 
of Wien and Rayleigh in the limiting cases of short and long wavelengths re- 
spectively, while their experimental measurements favoured it over all others. 

However, the derivation was largely phenomenological since the entropy had 
not been calculated from first principles. Therefore Planck (1900Vb) next 
sought to follow the "Boltzmann prescription" whereby entropy involves dis- 
order through the relation S = k In W, where W is the number of ways or states 
in which the physical system may be realized. This use of this famous relation, 
never given by Boltzmann himself yet inscribed on his tombstone, is the first 
occurrence known to the writer. Since the total radiation energy is distributed 
among all the radiation resonators involving a number of frequencies, it is first 
necessary to consider the distribution of an energy E among the N resonators 
with the same frequency v. Now if E were infinitely divisible the distribution 
would be possible in infinitely many ways. Hence Planck considered - and this is 
the essential point of the whole calculation - E to be made up of a fixed number 
P of finite equal parts, or energy elements, each of magnitude e = hv, where h 
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is a new constant of nature, so that E = Pe = Nw, where w is the average en- 
ergy per resonator. This therefore represents the first appearance of the Planck 
constant, both as to its symbol and definition. Planck also proceeded to use 
combination theory to calculate the total number of possible complexions : 

W = {N + P - 1)1 /(N - 1)!/P! » (N + P) N+P /N N /P P 

on making use of the Stirling formula whereby for large n, n! n n . He did not 
however give an expression for the entropy S, but gave his final result for the 
spectral distribution in its definitive form : 

u{v,T) = (87r^ 3 /c 3 )/[exp(^/fcT) - 1] 

where now the Boltzmann constant k appears explicitly as well as the new 
(Planck) constant h. From the formula, using the measurements of Kurlbaum 
(1898) and of Lummer and Pringshcim (1900), the values of h and k were 
determined to be h = 6.55xl0 -27 erg sec and k = 1.346xl0~ 16 erg grad. 

Very shortly afterwords, Planck (1901a) completed his derivation with details 
of the calculation of the entropy giving : 

\nW = N[{1 + P/N) ln(l + P/N) - (P/N) ln(P/N)] 

so that, since P/N = w/e, the entropy per resonator becomes : 

s = S/N = fc[(l + w/e) ln(l + w/e) - (w/e) \n(w/e)} = f(w/e) 

which for small values of w/e (small w and large e or v) leads to the earlier 
(Wien) result. Now it follows that : 

As i 'Aw = 1/T = (k/e) ln(l + e/w) 

giving : 

w = e/{exp(e/kT) - 1} 

Then since : 

u(v, T) = 8tt(^ 2 /c 3 )w 
and, according to the Wien Displacement Law : 

u{v,T) = v 3 f v {v/T) 

where f v is a still undetermined function : 

w = v[(c 3 /8ir)Mv/T)} 

Thus : 

1/T = As/Aw = v- 1 {8n/c?)f- 1 {w/v) 
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so that, on integrating : 

s = F(w/v) 

and therefore, in order to agree with the result that s = f(w/e), it follows that 
e = hv giving : 

w = hv/[exp(hv/kT) - 1] = (hc/X)/[exp(hc/kTX) - 1] 

where h is again the Planck constant. Thus one finally obtains the Planck 
formula : 

u(v,T) = {8irhv 3 /c 3 )/[cxp{hiv/kT) - 1] 

u(A,T) = (8irhc/\ 5 )/[exp(hc/kT\) - 1] 

The integrated radiation energy density is now u(T) = aT 4 with (Planck, 
1901a) : 

a = (8irk 4 /h 3 c 3 ) / z 3 /(exp z - l)dz = 8ir 5 fc 4 / '15c 3 h 3 
Jo 

in which the integral evaluates to 3!£(4) = 7r 4 /15 giving a = 7.565xl0~ 16 m~ 3 i^~ 4 
From the expression for w, giving s = fc[z/(expz — 1) + ln{exp z/(exp z — 1)}], 
where z = hv/kT, the integrated entropy density is S(T) = (4/3)aT 3 F, as has 
already been derived from thermodynamics. 

The maximum of the wavelength distribution occurs for x m — hc/kTX m 
where : 

x m = 5[1 - exp(-x m )} = 4.965114... 

while the maximum of the frequency distribution occurs for y m = hv m /kT 
where : 

y m = 3[1 - exp(-i/ m )] - 2.821439... 

giving v m \ m jc — y m /x m = 0.568252..., indicating that A m < c/v m and v m < 
c/X m . Hence also : 

u m /T = {k/h)x rn = 5.8785xlO w HzK- 1 

X m T = (hc/k)/y m = 2.8§77§xlO- 3 mK 

In the low frequency/long wavelength limit hv/kT — hc/kTX <C 1 the Raylcigh 
formula is obtained, on expanding the exponential in the denominator : 

u{v,T) = {8nhv 2 / c 3 )kT 

u(X,T) = 87rfcT/A 4 

while in the high frequency /short wavelength limit hv/kT — hc/kTX 3> 1 the 
Wien formula is obtained, on ignoring the 1 in the denominator : 

u(v,T) = (8tt/i^ 3 /c 3 ) exp(- hv/kT) 
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u(X,T) = (8tt/ic/A 5 ) exp(-hc/kTX) 



Further experimental measurements by Lummer and Pringshcim (1901) showed 
that the values of X m T were in excellent agreement with the predictions of 
Planck. Finally, Planck (1901b) returned to the question of the entropy of radi- 
ation, to demonstrate that the expression obtained from his theory satisfied the 
maximization principle. Thus in his stumbling search for the blackbody radia- 
tion spectral distribution, Planck, in forging the hitherto missing link between 
radiation and thermodynamics, not only closed the era of classical physics but 
also opened up the era of quantum physics. 

5 Biographical Note 

Max (Karl Ernst Ludwig) Planck was born on 23 April 1858 in Kiel, and died 
on 4 October 1947 in Gottingcn. He attended the universities of Munich and 
Berlin (in the latter he studied under Kirchhoff and Helmholtz), and received 
his doctorate from the University of Munich in 1879. Subsequently he taught 
in the universities of Munich (1880-1885), Kiel (1885-1889) and Berlin (1889- 
1928). In 1894 he became a member of the Prussian Academy of Science, and 
the permanent secretary of its Section of Mathematics and Natural Science in 
1912. He was awarded the Nobel Prize for Physics in 1918. In 1930 he became 
president of the Kaiser Wilhclm Society in Berlin, later to be known as the Max 
Planck Society. 
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